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We clarify the domain needed for the mixing angles in 
three flavor neutrino oscillations. By comparing the ranges of 
the transition probabilities as functions of the domains of the 
mixing angles, we show that it is necessary and sufficient to 
let all mixing angles be in [0,7r/2]. This holds irrespectively 
of any assumptions on the neutrino mass squared differences. 
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I. INTRODUCTION 

With the growing realization that neutrinos might 
have non-zero masses, studies of neutrino physics data 
in terms of neutrino oscillations, with three (or more) 
flavors, are becoming more and more popular. 

The easiest way to describe the neutrino oscillations is 
to use a set of two flavor oscillation schemes. However, 
since we observe three neutrino flavors in nature, a more 
natural approach is to use a three flavor oscillation model, 
and there are reports of different results obtained, when 
using two or three flavors [Qj|] . 

Much work has been done on neutrino oscillations us- 
ing two flavor models. When using a three flavor os- 
cillation model, one simplification often used is to make 
assumptions of the mass squared differences to enable the 
problem to be treated using much of the results coming 
from the use of two flavor oscillation models. 

In a two flavor oscillation model it is natural to let 
the mixing angles be in [0,7r/4]. By assuming the mass 
squared differences to be in certain ranges that allow the 
three flavor oscillation model to be simplified, using re- 
sults from the two flavor oscillation model, this domain 
of the mixing angles is often inherited to the three flavor 
oscillation models. Several authors that have made such 
assumptions about the mass squared differences have re- 
stricted the mixing angles to [0, 7r/4] Q. Other authors, 
however, have found solutions with mixing angles larger 
than 7r/4 when using a three flavor oscillation model M . 

To clarify the situation we here perform a detailed in- 
vestigation of the domain needed for the vacuum mixing 
angles without any speciflc assumptions and show that 
for a correct analysis to be possible, it is necessary to let 



the domain be < 6*^ < n/2 for aU i = 1, 2,3 g. This 
holds even if assumptions on the mass squared differences 
are made. 

Our paper is organized as follows. In Section 2 we de- 
scribe our formalism. In Section 3 we analyze the domain 
needed for three flavor oscillations, and in Section 4 we 
discuss the results. 



II. FORMALISM 

For the present analysis we will use the plane wave 
approximation when describing the neutrino propaga- 
tion. We will describe the neutrino flavor states pro- 
duced by charge current weak interactions as a mixing of 
neutrino mass eigenstates. For simplicity we neglect any 
CP-changing phases. We then write the flavor states 
i^ct), a = e, fi,T, as a linear combination of the mass 
eigenstates [vi), i = 1,2,3, 



|i^a) =^^*,k»), a = e,/x,T. 



(1) 



Here U is a unitary mixing matrix, given by the stan- 
dard representation of the Cabibbo-Kobayashi-Maskawa 
mixing matrix as 



C2C3 
U = I — S'aCi — S1S2CZ 
Si S3 — S2C1C3 



S3C2 S2 

C1C3 ~ S1S2S3 S1C2 

— SiCs — >S'2S'3Ci C1C2 



(2) 



with Ci = cos9i and Si = sin^i, i = 1,2,3. Since there 
are no CP phases, the matrix elements are real and we 
have C/*(j = Uaa for a = e, /z, r and a = 1, 2, 3. 

The transition probability for a neutrino changing from 
i^a to vp is given by 



Pal3 = S, 
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Here Am? 



I?, E is the 
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where a, P = e,ii,T. ^^^^^ —■^ij — •■-i ■■-j' 
relativistic energy of the neutrinos and L is the source- 
detector distance. 

The mass squared differences are not all independent 
in a three neutrino model, but fulfill by definition the 
relation 



Am2i + Am 



Amj^ 



0. 



(4) 
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From Eq. (||) it follows that Paf} = Ppa- 

In what follows we shall, without loss of generality, as- 
sume that the masses are ordered so that mi < m,2 < Ti3, 
and the mass squared differences Attj?- are non-negative 
when i > j. 

From the unitarity condition we have three equations 
that the transition probabilities must fulfill: 






(5) 
(6) 

(7) 



Since we want to investigate the mixing angles 0^, we will 
for convenience introduce the auxiliary 'angles' 



AtooiL , Am%nL 

$1 = — -^f— and $2 - 



4S 
Using Eq. (0), this gives 



4^ 



4^ 



(8) 



(9) 



The auxiliary angles are both positive by definition. To 
introduce the auxiliary angles makes sense, because we 
do not here care about the mass squared differences. The 
values of $1 and $2 vary with L/E, which is something 
that is determined by the experimental setup. The mass 
squared differences only determine how fast the auxiliary 
angles vary. 

We can now write Pap as 



P. 



a/3 — OafJ 



4[/^^sin2$i-H/2^sin2<D2 



+ /3^sin2(ci>i + $2)], 
where /^^ = f^piOi, 02,83) is given by 

faf3 — UalUpiUa2Up2, 
fafS — Ua2Up2Ua3Up3, 
faP — UalU/3iUa3U/33. 



(10) 



(11) 



III. ANALYSIS 

To determine the necessary and sufHcient domain for 
di, i = 1,2,3, we first show that it is sufficient to take 
6i G [0,7r/2], and then continue to look whether it is 
possible to restrict the domain further. 

We first notice, that due to unitary, it is sufhcient to 
look at the diagonal probabilities, P^e, Pf^fi, and P^-^. 
The nondiagonal probabilities Pe^ , Per , and P^t , can be 
calculated from the diagonal ones using Eqs. (p[)-(R). 

The diagonal probabilities are given by 



Paa = 1 - 4 ^ ^ C/2, [/2^. siu^ "^'^ ^ 



i=l j=l 
i<3 
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a = e, fj,, T. 



(12) 



Since the UaiUaj appear everywhere as U^^U^j, every 
sin 9i and cos 9i will appear as squared, causing the signs 
to disappear. It is therefore always possible to choose 
an angle between and tt/2 to parameterize the mixing 
angles. A sufficient domain for their variation is therefore 
given by [0,7r/2]. 

In order to see if it is possible to reduce the domain for 
di further, we want to find a way to handle the problem 
that the Pqq's are functions of five variables, two mass 
squared differences and three mixing angles. We will do 
that by fixing $1 and ^2- For every fixed $1 and $2 
we then have the PaaS as functions of the three mixing 
angles. This means that we have solved the problem of 
having too many parameters at the cost of getting a lot 
of functions instead. For each pair of $1 and $2 there is 
a function that depends on three parameters, the mixing 
angles. 

Next, it turns out to be practical not to explicitly look 
for the domain needed for the mixing angles, but instead 
check if there is any difference between different domains 
of the angles. If we look at one of the Pqq's, it will take 
some range of values for some fixed domain of values of 
the 0i's. If there are some values it cannot take when the 
angles are restricted to a smaller domain, the restriction 
makes a difference. We will therefore look at the range 
of Pqq, denoted R{Paa)- 

Considering one probability at a time and without any 
restrictions on the parameters, the range of Pqq is of 
course P(Pqq) = [0,1]. However, for fixed values of $1 
and $2, the range R{Paa) can be smaller. To see this, 
take L — 0, giving $1 = $2 = 0, as an example. This 
leads to Paa — 1 for all possible fl^'s, which is obviously 
a smaller range. 

To simplify the forthcoming discussion we introduce 



^L=4[/^„sin2$i+/2^sin2$2 
+ /Lsin2(ci>i+$2)]. 



(13) 



We do this because all f^^, fc = 1, 2, 3, have at least one 
factor sinOi for some z, which means that 



P' = when 61 = 62 = 93 = 0, 



(14) 



and we thus obtain a range of the type [0, a;], where a; < 1, 
for all Pqq's. We see that the interesting value is the 
upper limit of P(Pqq). 

From now on we will always consider $1 and $2 as 
fixed when calculating R{P^^), making R{Paa) a func- 
tion of $1 and $2- This means that we take fixed $1 and 
$2 and then calculate the range of the function, which 
now has the three mixing angles as parameters. Each 
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FIG. 1. Pee, *i = 0.7 and $2 = 2.0. Note that P^^ does 
not depend on 611. 



choice of $1 and $2 will give a specific range. The up- 
per limit of the range can obviously change for different 
domains of the 6'i's. Introduce the notation 

RiPDabc = { PU0i,02,e3) : di e [0,a], 

92e[O,b],03e[O,c]} (15) 



We now calculate R{Paa)abc, abc 
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the /^ s are 'nice' functions, that vary slowly, this is not 
a problem to do numerically. Fig. m shows an example 
of how smoothly P^g varies with 9i. 

The result is, that for some values of $1 and $2 we 
have 



i?(pL)fff y^^(^L)^ 



(16) 



This means that for some values of $1 and $2, there 
are values of the functions P^„ that can be obtained by 
letting 0i e [0, 7r/2], but not be when 0i G [0, 7r/4]. 
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FIG. 3. Plot of R{P^^)^^^ - R{P^, 



In the same way, one can check what happens when 
only one or two mixing angles are restricted to [0, 7r/4]. 
The result is that it is possible to restrict either di or 62, 
but not both, to the smaller domain without getting any 
differences in the ranges of the Pea's. 

To visualize this result we define the difference of 
R{Pabc) ^^'i ^{P'a'b'c') ^^ ^^^ difference of the upper hmit 
of the ranges 



R{P'aa)abc " R{P'aa)a'~ 

= max (Paa) - 

yi,2,3<O,0,C 



inax /P^J. (17) 

'l,2,3<a ,0 ,c' 



We can now make a contour plot of this difference with 
$1 and $2 on the axes. The result is displayed in Figs.g- 
Q. In the empty areas of the Figs, there is no difference 
between the domains, but inside each curve there is a 
difference that gets larger as we move inwards. The plots 
have a two-fold rotational symmetry around the point 
$;^ = $2 = 7r/2. 

For example, in R{P'^^) we see from Fig. 0, that the 
maximum difference occurs when $1 ~ 0.7 and $2 — 2, 
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FIG. 2. Plot of R{Pi 



RiPDf. 




FIG. 4. Plot of R{Pr^ 



(or <I>i ~ 2.4 and $2 — l-l)- 

If we look at the probability at that point, we have 

P^, = 4 [ {C2C^S3C2f sin^ $1 + (53^2^2)' sin^ $2 

+ (C2C352)2sin2($i + $2)]. (18) 

Since this equation only depends on 62 and ^3, its range 
can easily be shown in the plot displayed in Fig. H. It is 
clear that the range is larger when allowing ^3 > 7r/4. 

To see how large the difference is, we have given the 
maximum difference for each of the P'^^ in Table |. It can 
readily be seen that the maximum differences are large 
enough to be of importance. 



i'^aa 


Maximum difference 


Pee 


^0.19 


"mm 


^0.05 


Prr 


Ri0.19 



TABLE I. Maximum difference between R{P^a) ^ ^ ^ and 

-R(-Pac)f f f • 

For the nondiagonal probabilities the differences are 
even larger. 



IV. DISCUSSION AND CONCLUSION 

To see which areas of the plots that are relevant, we 
look at the definition of $1 and $2 in Eq. (||). The 
two mass squared differences have fixed physical values 
(although largely unknown at present). The quantity 
L/E can be used as a parameter for $1 and <I'2- Varying 
L/E corresponds to a line in the plots in Figs. 0-0 modulo 
w. An example of how the line may look like is given in 
Fig. I 

The slope of the line depends on the ratio between the 
two mass squared differences and the size of the area it 
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covers depends on what values L/E can take. We see 
that the slope of the line is of little importance when we 
check if it intersects a marked area in Figs. |2|-^. This 
means that, even when making assumptions of the mass 
squared differences to simplify the analysis, it is necessary 
to allow the mixing angles to be in [0,7r/2]. 

When considering a line as in Fig. ^, one should re- 
member that the length of the line depends on the range 
of values that L/E can take in an experiment. In a par- 
ticular experiment this may be a very short linesegment, 
but this does not imply that we only have to consider 
this short line. By making the experiment better, or us- 
ing other experiments measuring the same probability, 
the line will be longer and thus any imaginable extension 
of the experiment should be taken into account when con- 
sidering the line in Fig. 0. This means that in practice 
the line will always cover the whole plot. 

The analysis above showed that it is possible to di- 
minish the domain of either 9i or 62 without getting any 
difference of the ranges of the PaaS. It is, however, not 
at all clear which of them one should choose to diminish. 
As different choices can lead to different solutions, both 
of them should be left to be in the domain [0, 7r/2]. 

It is possible to show, using the same technique as 
here, that as soon as one tries to diminish the domain 
of variation of all three mixing angles from [0,7r/2], it 
makes a difference in the ranges. 

In conclusion, when using a model for three flavor os- 
cillations as described briefly above, it is necessary and 
sufficient to allow all of the mixing angles to be in the do- 
main [0, 7r/2] when fitting experimental data, even when 
making assumptions on the mass squared differences. 
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